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Abstract 

Multiport generalizations of beam splitters are the key component in multipath interferometers, 
which are important in a range of quantum state engineering and precision measurement schemes. 
Here we propose a straightforward method for implementing multiport devices for atoms trapped 
in optical ring lattices. These devices are interesting as atoms have certain properties (such as 
mass) that photons do not and the ring configuration makes them useful for applications such as 
precision gyroscopes. We discuss how they could be employed in useful measurement schemes and 
investigate how practical considerations limit the size of the devices that can be achieved by this 
method. 

PACS numbers: 03.75.-b, 03.75.Lm, 03.67.-a 
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I. INTRODUCTION 



Multiport generalizations of beam splitters enable the manipulation of quantum states 
which may be important for quantum information processing on networks. They are also 
a key element in multipath interferometers and so have great potential in a range of mea- 
surement technologies. Considerable progress has been made over the past decade towards 
implementing and understanding these devices. Theoretical work has devised schemes to 
create multiport devices and analyze some of their useful properties js,!^!. Experimental 
work using photons has demonstrated the operation of beam splitters with three inputs and 
outputs ('tritters') as well as four inputs and outputs ^quarters') by clever arrangements 



of ordinary two-port beam splitters and phase plates 



J]. Similar devices have also been 



proposed ^] and implemented in systems of bundled optical fibres. 

These devices become more powerful as the number of input and output ports is in- 
creased. However, for larger systems, the experimental configuration required rapidly be- 
comes labyrinthine. In general, for a device with 5* inputs and 5* outputs, S{S — l)/2 beam 
splitters and S{S + l)/2 phase shifts would be required, i.e. a total of 5*^ optical elements 
The complexity of implementation is perhaps even more problematic for atomic systems 
where beam splitting is a dynamic process typically involving the raising and lowering of 
potential barriers or the careful application of Bragg pulses. In this case, a multiport beam 
splitter would involve a complex sequence of operations and the experimenter would need to 
be able to address lattice sites individually - an issue that has caused considerable problems 
when trying to use optical lattices for quantum information processing. 

In this paper, we demonstrate a scheme for implementing atomic multiport devices that 
goes someway towards overcoming these problems. In particular, it requires no additional 
equipment or steps than a standard two-port beam splitter. We begin by developing the 
theoretical scheme for an atomic tritter and follow with its extension to a device with a 
general number of ports. An interesting asymmetry is found in the behaviour of devices 
with an odd or even number of input ports. We then highlight how this scheme could be 
used for measurements. Finally, we take into account various practical considerations and 
see that our scheme is likely to be limited to splitters with about five input (and output) 
ports. 



2 



II. THE SCHEME 



The physical system we consider consists of an optical lattice of 5* sites in a 'ring' con- 
figuration with atoms trapped at the potential minima. This can be described by the 
Bose-Hubbard model with Hamiltonian, 

5-1 5-1 5-1 

j=0 j=0 j=0 

where aj is the annihilation operator for an atom at site j and the ring geometry means that 
O'S = (^0- The parameters J and V are the coupling and interaction strengths respectively, 
and ej accounts for the energy offset of site j. In general, we take the zero point energy to 
be the same for each site and so can ignore these energy offset terms. However, making one 
or more of the energy offsets non-zero for a fixed period of time can be a convenient way 
of imprinting phases on individual sites [7|. Initially we take the potential barriers between 
the sites to be sufficiently large that we can ignore tunnelling. This now is the starting 
configuration and each lattice site corresponds to an input port of our multiport device. 

The first step is to rapidly reduce the potential barriers between the sites in such a 
way that the sites still remain separate but are strongly coupled due to tunnelling. We 
want to do this rapidly with respect to the tunnelling time, but slowly with respect to the 
energies associated with excited states to ensure our system remains in the ground state. 
The adiabaticity criterion that ensures no excited states are populated is given, in general, 
by the phonon excitation spectrum {ujk}. This can be found using Bogoliubov theory and 
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has the well-known form 



uJk = \ 4:Jsm 
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where N is the total number of atoms, 5* is the number of lattice sites, and the index k runs 
over values k = 0,1, 2, {S — 1). Experiments have already successfully demonstrated this 



io| 



separation of timescales by ramping the optical intensity on a timescale of about 20ms 
The limitations imposed by this condition are discussed further in Section |Vl 

The coupling between wells now dominates over the interactions and the Hamiltonian ([T]) 
can be written in the simple form, 

5-1 
j=0 



Of course this makes the approximation that the interaction energy will be negligible 
compared with the coupling energy. In practice, the interactions will not be strictly zero. 
However, they can be made very small with respect to J by, for example, making use of 
Feshbach resonances to tune the scattering lengths llj. For now, it is helpful to consider 
the case where we can ignore V . We consider the effects of interactions in Section IVl 



A. Example: The Tritter 



We begin by considering the case of three lattice sites in a ring (see Figure [T]). A similar 
configuration has been achieved experimentally by trapping atomic Bose-Einstein conden- 
sates (BECs) in the optical potential created by the diffraction of a laser beam by a liquid 
crystal spatial light modulator [12]. This modulator allows arbitrary three-dimensional trap- 
ping potentials to be achieved, which have the added advantage of being able to be varied 
smoothly with time. Another promising possibility for creating the ring potential, is to 
interfere a L aguerre- G aussian (LG) laser beam with a plane wave co-propagating along the 
2;-direction ISj. By retro- refiecting this combined beam, a standing wave can be formed that 
consists of a stacked array of disk shaped traps along the z-direction. By controlling the 
tunnelling between the disks and making it much smaller than the corresponding tunnelling 
within each ring, one can implement an array of effective 1-D ring lattices. In both these 
cases, the rate of tunnelling between the sites in a ring can be controlled simply by adjusting 
the intensity of the trapping laser light. 

The Hamiltonian ([3]) describing this system is. 



H = — J ^ajfli + a[a2 + a\ao^ + h.c. 



(4) 



This can be diagonalized in the basis: 
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If the system is allowed to evolve for time, t, the ao mode acquires a phase of —2Jt, while 
the «! and a2 modes each acquire a phase of Jt. If the barriers are then raised on a similar 



timescale to their lowering (i.e. quickly with respect to the tunnelling time, but slowly with 
respect to the energies associated with excited states), the atoms are 'frozen' in the lattice 
sites oq, ai, and 02. The overall output operators {^0,^1,^2} are given by: 
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, and for later reference, we have defined the 



where f2i = e + 2e~*'^ , f22 = e'""" — e 
overall operation as R^. We see that the output modes are identical to the input modes 
when t = 0, as we would expect. For a balanced fritter, we need the output modes to be an 
equal superposition of the input modes, i.e. = \^l2\- This occurs when t = 2tt/9J, for 
which value we have (ignoring any irrelevant overall phase). 
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To summarise, the steps in implementing a balanced fritter are: 

1. Rapidly reduce the potential barriers separating the lattice sites. 

2. Allow the system to evolve for time t = 27r/9J. 

3. Rapidly raise the potential barriers. 

Each output is an equally-weighted superposition of all the inputs. One may notice, 
however, that the phases between terms in these outputs are different from those com- 
monly quoted in the literature (e.g. in [2], Q]). For many purposes, this does not mat- 
ter. For example, we can take the output modes of a 50:50 beam splitter to be either 
{(ao + ai)/-\/2, {ao — ai)^/2} or {{ao + iai)/\/2, (iaQ + ai)\/2} without fundamentally chang- 
ing the results. If, however, we do require a particular form of the phases for the output 
modes, this can always be achieved by imprinting phases on individual lattice sites before 
and after the lowering and raising of the barrier. In practice this could be achieved by 
applying energy offsets, ej, to the lattice sites for some fixed time, as discussed above. 

This procedure not only achieves the goal of a multiport beam splitter for atoms, but 
is much simpler than schemes that combine phase shifts with a complicated network of 



beam splitters. Our scheme only requires a lowering and raising of the lattice potential and, 
importantly, requires no more operational effort than an ordinary beam splitter for atoms 
14| . This bodes well for the possibility of scaling the scheme up to larger systems. 



III. LARGER DEVICES 



We now extend this scheme to the general case of a lattice ring with an arbitrary number 
of sites, 5* (in Section |V] we will determine up to what values of S this scheme is valid exper- 
imentally). Following the same procedure as for the tritter we lower the potential barriers 
so that coupling between the wells dominates over the interactions and the Hamiltonian 
describing the system is given by ([3]). This can be written in the diagonalised basis as, 

5-1 
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The system is now allowed to evolve for time t during which the S modes each acquire a 
phase of —2 cos (^) Jt. The potential barriers are then quickly raised and the output is 
given by. 
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where the unitary matrix is given by U^j 

As before, the condition for a balanced multiport splitter is that each output is an equally 
weighted superposition of all the inputs, i.e. \VLi\ = \ fl2\ = • • • = \^s\- We now wish to see 



if our procedure can achieve this for a general S ignoring, for now, factors that may hmit 
its experimental feasibility. 

A. Producing a balanced splitter 

We can determine if a balanced splitter can be achieved for a device with S sites by 
plotting the elements of Rs as a function of Jt and checking whether there are times for 
which = \Q2\ = ■ . ■ = \^s\- Figure [2] shows these plots for 5 = 3,4,5 and 6 respectively 
in the time range t = —>■ 2-71/ J. We see that all the Q intersect at t = 27r/9J for S = 3 
and at t = n/AJ for 3 = 4, but there is no exact crossing for S = 5 and 6 in this range (or 
indeed for larger values of S). However, looking in the time range t = —>■ 5000/ J we find 
for systems with odd values of S, up to S* = 9, there is an intersection of all the Q to within 
an error of 1%. It is convenient to define a measure, Xy of how good the intersection is, as 

x^|:(|n.|--L)\ a,) 

This is simply the sum of the squares of the deviation of each value of from its ideal 
value (i.e l/v^) - as discussed above, we are not concerned by the phase of the f2 values. 
If, on average, each value of differs from by 1%, x has a value of 10~^. We will 

use this value (i.e. x ^ 10"'*) as a useful criterion for when a balanced splitter has been 
achieved. 

Using this criterion, we find balanced splitters for S = 5,7 and 9 at evolution times 
t ~ 5.2tt/J, 8871/ J and ITTvr/J respectively. This intersection for 5 = 5 is shown in 
Figure O The time required becomes longer because there are more values of Q to match. 
Interestingly, no intersections for which x ^ 10~^ are achieved for even values of S in this 
time range. We will discuss possible reasons for this in the next section. We should note 
that while we do not see intersections that satisfy x ^ 10~^ for even values of S in the time 
range used, this does not rule out the possibility at a later time. One exception is 5* = 6 
because all the values of Q in this case repeat with a period of 27r and, since there is no 
intersection during this period, there will never be one. The same may be true for some 
other, larger S. However, so long as x for a given 5* is not periodic in time, if we wait long 
enough the balanced splitting criterion should be able to be achieved. Of course, this may 
not always be experimentally expedient and could be longer than the coherence time of the 
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system. 

To summarise, devices with even and odd numbers of sites are potentially both capable 
of producing balanced multiport devices. However, the even cases usually require a much 
longer evolution time than the odd ones. We now investigate the reasons for this difference 
in behaviour. 



B. Evolution while the barriers are lowered 



For simplicity we will begin by considering the evolution of a state containing only a 
single particle. However, our results also apply to more general states. An initial single 
particle state can be written as, 

5-1 5-1 

|$) = ^A,at|0) = ^i?,aI|0), (13) 

j=0 k=0 

where Aj and Bk are the amplitudes in the two bases aj and respectively, and are related 
by a Fourier transform. Using fllOj) and fll3p we can write the evolution of the state 
while the barriers are lowered as, 

^ ^ ''27rA;\ i27rkj 



fc,J=0 



i2Jt cos 



S S 



a]\0). (14) 



Taking the case where the single atom is initially in site 0, i.e. Aq = 1 and A„ = for 
n = 1, 2, . . . , S* — 1, for an even value of S we can write, 

s-l 5/4-1 

\m) = E^"'^' E mk,s,t)a]\o) (is) 

j=0 k=l 

where A(j, k, S,t) is some real- valued function that depends on j, k, S and t. This can be 
written more simply as, 

5-1 

|$(i))=EA(j,^,t)e^^^/^a]|0), (16) 

j=0 

where A{j,S,t) is also a real-valued function. Importantly from (fT6|) we see that adjacent 
sites always have a phase difference of ±7r/2. This means that the rate at which atoms can 
flow between sites has certain fixed values since the flow velocity is given by, v = {h/m)'V(j) 
where V0 is the phase gradient or phase difference. 
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By contrast, for an odd number of sites the wave function after evolution for time t can 
be written as, 

S-l ( 5/4-1 \ 

= y\e^''^e'^'l^ A{j,k,S,t) 1 a]\0) 




+ e^''^/'Aij,S,t)J 410). (17) 

Here the phase at each lattice site varies continuously with time and so the phase differences 
between adjacent sites, and hence the velocity of flow, are not constrained in the same way 
as the even case. The additional constraint imposed by the symmetry of the even case on 
the allowed flow rates is one reason why it is more difficult for devices with an even number 
of sites to achieve a state where their outputs are in an equal superposition of their input 
modes 2^. We have taken the simple example of a single particle state, however the same 



argument holds for general input states: the only change is the form of the real-valued 
functions. 

Another difference between the cases of even and odd numbers of sites can be found 
by considering the direction of the flow of atoms around the ring. In the case of even S 
the amplitude of each site is given by e*'^-'/^A(j, S,t), and so the phase is always ±1 or ±i. 
The function A(j, S, t) is smoothly varying and so changes sign only when it passes through 
zero. Consequently the phase difference between two sites, and hence the direction of flow, 
changes only when the amplitude of one of the sites is zero. The same is not true for the 
devices with odd 5* as their phases change continuously with time. This allows the direction 
of flow between the sites to change smoothly with time. This difference in the evolution of 
systems with an even number of sites and those with an odd number of sites again suggests 
it will be more difficult for those with even S to achieve the required balanced output. 

The final difference is seen when we look at the number of distinct elements required 
to compose the operators, Rs- Instead of the number of different Qs equalling the number 
of sites, S, there are {S + l)/2 different Qs required for odd S and {S + 2)/2 different fls 
for even S. This is due to the symmetry of the system and it means for even S there are 
more values of Q to match to make a balanced splitter. Again, this difference supports our 
observation that it is more difficult for systems with an even number of sites to produce a 
balanced splitter. 
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IV. INVERSE TRANSFORMS AND UNBALANCED SPLITTERS 



Now that we have understood the basic operation of multiport sphtters, we would hke 
to consider how they can be combined to create useful quantum devices. One simple, 
yet important, such device is an interferometer. In general, an interferometer consists of 
a beam splitter, a phase shift, and then an inverse beam splitter. In standard two-path 
interferometry, the inverse beam splitter is replaced with a normal beam splitter. The 
reason this works is that a 50:50 beam splitter can be thought of as a vNOT operation. 
We can see this because a Mach-Zehnder interferometer with no phase difference between 
the two paths (i.e. two 50:50 beam splitters in succession) gives an output state that is the 
same as the input but with the ports swapped, i.e. a NOT operation. So two beam splitters 
in succession will return us to the original state so long as we make a trivial swap of the 
labels of the output ports. For multipath interferometers, however, the same is not true: we 
cannot simply replace the inverse splitter with an ordinary splitter and relabel the modes. 
So, to be able to implement a multipath interferometer, we need to be able to implement 
inverse multiport devices. 

It turns out that three successive applications of a tritter, is equivalent to the identity 
and leaves the original input state unchanged. This can easily be verified using Equation 
(IH]). This means that the inverse operation of a tritter, R^, is simply or, equivalently, 
a tritter operation where the state is allowed to evolve for twice as long with the barriers 
lowered (i.e. t = 47r/9J rather than t = 27r/9J). The inverse operation can, therefore, be 
implemented just as easily as the original transform. 

(S—l) 

More generally, the inverse operation of a device with S sites is simply Rg . We have 
numerically checked this for S = 4,5, 7, and 9 [23] by calculating how close RsR^g~^^ = Rg 
is to the identity operator. A useful measure of this is simply to sum up the modulus squares 
of the leading diagonal of Rg and divide by S since this effectively gives the fidelity of the 
output state relative to the input (averaged over all input ports). The values we found for 
this measure for S" = 4, 5, 7 and 9 were 1, 0.96, 0.91 and 0.75 respectively. We see that the 
inverse splitter is degraded as 5* is increased. This is not surprising because, if our timing 
does not give an exact balanced splitter, when we multiply this time by {S — 1) to give the 
time of the inverse splitter the imperfection is multiplied. 

So far, we have only discussed splitters that have balanced outputs and are therefore the 
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multiport generalizations of 50:50 beam splitters. For many applications we may not want 
the outputs to be balanced. For example, we may want a device that coherently skims off 
only a small fraction of an input state and redistributes it between the output modes. Our 
scheme holds great potential for producing such devices simply by changing the value of 
the evolution time. Since the coherent amplitude in each output mode depends on t, this 
allows us to obtain multiport splitters with different 'reflectivities'. All this can be achieved 
without changing the experimental set-up - only the timings of the steps. It is important 
to note, however, that not all unbalanced outputs can be achieved by this method. For 
example, in Figure [2] we see that not all the values of Q are independent. 

We now consider how these multiport splitters can be combined into useful devices. A 
simple, but illuminating, example is that of a three-path interferometer. Suppose we have 
three lattice sites in a ring configuration and start with atoms in one site, oq. The 
interferometer is implemented by a fritter operation, R^, a phase shift on the lattice sites, 
and then an inverse fritter, R^^. If we take the case that there is a linearly varying phase 
shift of (p between adjacent lattice sites the overall transformation of the interferometer is. 



(18) 



The mean number of atoms at each final lattice site (normalized by the total number) is 
then given by: 
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3 + 4cos(0) + 2cos(20) 
3 + 4 cos(0 - 27r/3) + 2 cos(20 + 27r/3) 
\^ 3 + 4 cos(0 + 27r/3) + 2 cos(20 - 27r/3) J 

These mean numbers of atoms are plotted as a function of (p in Figure HI As we might expect 
from such an interferometer, a measurement of the final population in each of the three lattice 
sites allows us to uniquely determine the value of (p (modulo 2tt). Importantly, we note that 
this three-path interferometer requires no more effort than a two-path scheme for atoms (and 
indeed a scheme with many paths would also require no additional equipment or steps). This 
is in contrast with existing schemes jl| that combine beam splitters in complicated arrays 
and illustrates the power of this scheme. 
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V. PRACTICAL LIMITATIONS TO THE NUMBER OF PORTS 



So far in discussing the general scheme for implementing an atomic multiport device we 
have neglected a number of important physical processes. This has allowed us to treat an 
idealised case and has been useful for understanding its key features. However, in reality, 
these processes must be accounted for and will limit the applicability of the scheme. In this 
Section we highlight the key practical limitations to this scheme and assess their impact. 



A. Condensate lifetime 



One potential limiting factor is how long the system must be allowed to evolve to produce 
the balanced splitters since this may be longer than the lifetime of the condensate. In Section 
mil we showed that the evolution times increased with S. Figure [5] shows the value of Jt 
required to produce balanced splitters for 5* = 3, 4, 5, 7 and 9 and we see that in this range 
Jt scales as approximately the seventh power of S. This is quite a prohibitive scaling and 
in order to determine what values of S are feasible, we need to compare the evolution time 
with the lifetime of the condensate. 



Condensate lifetimes of about 10s have been measured [15|] and oscillation periods of 40ms 



161 ] ■ giving J ~ 80s ^ 



have been observed between wells in bosonic Josephson junctions 
Together these give an estimate of a typical maximum value of Jt ~ 800. For S = 3, 4, 5, 7 
and 9 the largest value of Jt is 1777r which is comfortably below 800. Therefore, the times 
required to produce balanced splitters with up to nine sites should all be experimentally 
achievable. Using the scaling of Jt with S above, for 5* = 11, Jt > 800 suggesting S* = 9 is 
the maximum number of sites the evolution times allow. 



B. Intensity fluctuations 

Since the time for which the barriers must be lowered depends on the tunnelling rate, 
any uncertainty in J will introduce imperfections into the scheme. This may concern some 
readers since it is well known that J depends exponentially on the intensity of the trapping 
light which is subject to fluctuations. However, just because J depends exponentially on 
the intensity does not mean that small fluctuations in the intensity will result in large 
fluctuations in the tunnelling rate. Indeed the exponential function varies only linearly for 
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exponents close to zero, and this is the regime in which this scheme operates. We can check 



this using an approximation for J found in 

Er 



hJ 



exp 





(20) 



where Er is the atomic recoil energy, Vq is the potential depth and h = 1 m [12] • To ensure 
we are in the strong coupling regime we take Vq = 2£'/j, a configuration that has been 



achieved experimentally 18|. For a small fractional fluctutation, 5 in the exponent due, for 



example, to intensity fluctuations in the trapping laser, this can be rewritten as. 



J ~ exp 
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Intensity fluctuations can be stabilised to around 0.1% p^], i.e. 6 = 0.001, which gives 
J ~ 0.9965J. So the uncertainty in the tunnelling rate due to intensity fluctuations is 
modest and unlikely to be the major limitation to this scheme. 



C. Interactions 



Up until now we have neglected the interactions between atoms by setting = 0. We 
shall now consider how non-zero interactions limit our scheme. Taking them into account 
the Hamiltonian describing the system is, 

5-1 5-1 

H = -jJ2 + + ^ E (21) 

j=0 j=0 

First, we shall consider the effect of these interactions on our fritter for different numbers 
of atoms. We will then determine how this effect scales with S. 

We measure the effect of interactions on the fritter by comparing the outputs of an inter- 
ferometer, composed of the fritter and the inverse fritter, with and without the interactions 
taken into account. For the purposes of our numerical calculations, we will consider the 
speciflc initial state consisting of atoms in site 0, i.e. |A^, 0,0). Figure [6] shows how the 
fldelity between the output states with and without the interactions varies as a function 
of interaction strength and for different numbers of particles, A^. As we would expect, the 
fldelity is unity when there are no interactions regardless of the value of A^ and the effects 
of the interactions become more signiflcant as A^ is increased. 
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To determine the scaling of the interactions with we find the value of VN/ J required 
to give a fidelity of 0.95 (when compared with the case without interactions) for different 
values of A^. From this we determine the relationship between the number of input atoms 
and VN/ J at this critical fidelity to be, 

'VN' 



J 



^0.85Ar-""7 (22) 

F=0.95 



which is shown in Figure [71 This relationship means the effect of the interactions, V/ J, on 
the system scales with as 0.85A^^^ °''. This is an approximately linear scaling and so the 
effects of interactions in the system are not too destructive. It means we require VN/ J to 
be of the order of 0.85 for F = 0.95. Therefore, if we were to input 10^ atoms we would 
need V^/J ~ 10~^. We can make V very small compared to J by using Feshbach resonances 
to tune the scattering lengths. It is possible to achieve scattering lengths smaller than oq 



for some BECs 



20l | meaning we can tune V/ J to be of the order of 10 ^, which is promising 



for the feasibility of the scheme. 

Things are less promising, however, when we consider the effect of interactions in larger 
devices. Determining V/ J at the critical fidelity of 0.95 for a fixed number of atoms (A^ = 5) 
for 5 = 3,4,5,7 and 9 allows us to see how interaction effects scale in this range of 5*. The 
results are shown in Figure [HI We see that the V/ J that still allows us to achieve a fidelity 
of at least 0.95 decreases rapidly with 5*. For 5* = 7, we require V/J^ 10""^, which is at 
the limit of what is experimentally achieveable. Interactions therefore appear to limit the 
current applicability of our scheme to about five sites. 



D. Timing errors 

So far we have made the simplifying assumption that we can exactly measure the time the 
barriers are lowered for. We will now investigate how sensitive the scheme is to inaccuracies 
in this timing for different numbers of atoms, A^, and different numbers of sites, S. We first 
consider a tritter (i.e. S* = 3). To determine the effect of an absolute time error, e, for 
different A^ we calculate the fidelity between the outputs of the device when Jt = r and 
when Jt = T + e, where r is the Jt required to produce the balanced splitter quoted in 
Section [ITI| for the particular input state |A^, 0,0). As with the interactions, we determine 
how the fidelity scales with A^ by measuring e at a critical fidelity of 0.95 for different values 
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of A^. These results are shown in Figure [9] and we see that the relationship between e at the 
critical fidelity and the number of atoms is given by, 

6^=0.95 ~ 0.16iV-o-50_ (23) 

The maximum time error that still allows us to achieve a fidelity of at least 0.95 decreases 
relatively slowly as increases. The same scaling of the absolute time error with was 
also numerically found for = 4, 5, 7 and 9. 

However, in an experiment it is the fractional error that we must consider since clocks 
generate fractional errors. We have determined these fractional errors for 5 = 3, 4, 5, 7 and 
9, when we input five atoms, by dividing e at F = 0.95 by r. These results are shown in 
Figure [TOl As expected the fractional error, e/r, that can be tolerated in these systems 
decreases with S. For S* = 9, e/r is of the order of 10~^. Accuracies of one part in 10'^ 
are easily achievable with many readily available pulse pattern generators suggesting that 
fractional time errors will not be the main limiting factor in this scheme. 

Our scheme relies on us being in the strong coupling regime where the tunnelling rate is 
larger than the oscillation rate. To ensure we satisfy this condition we take Vq = 2Eji and 
using fl20l) we find hJ ^ O.OGEr. Substituting values for ^^Rb and using A ~ lOOOnm we get 
J ^ 60Hz. We want to see how accurately we must be able to control the time given we are 
operating in a regime with this J. From Figure [3 we know the maximum absolute time error 
our system can tolerate is e ~ 0.07 for A^ = 5 for all 5* where e = J At and At is the time 
error in seconds. This means we require that we can control t to at least a precision given 
by At = 0.07/ J ~ 1ms. Therefore, we need a refresh frequency of about a few kilohertz 
which can be achieved using a ferroelectric liquid crystal spatial light modulator [21] and so 
we conclude that operating in the strong coupling regime does not limit our scheme. 

E. Rate of switching the intensity 

Another potential limit on S comes about from the rate at which the barriers are raised 
and lowered. As discussed above, we want this to happen quickly with respect to the rate of 
tunnelling between the sites, J, but slowly with respect to the phonon excitation spectrum 
in each trap. This ensures that the system remains in the ground state of each trap, and any 
imperfection in this switching will reduce the fidelity of the output state. This separation 
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of timescales has been experimentally demonstrated and is given by the ratio uJk/J., which 
can be calculated using Equation ([2]). In particular, we require uJk/J » 1 to ensure the 
scheme is successful. Putting different values of S into ([2]) and taking VN/ J to be 0.85 we 
find uJk/J> 1 for all values of 5 up to S* = 12. This gives a limit on S of 12. 



F. Loss of particles 

The final limiting factor we consider is spontaneous emission which we show to have little 
effect on the system. In fact, when we ignore interactions, the loss of a particle during the 
splitting process has no effect on the remaining particles. When particles are inputted into 
a fritter and one particle is lost halfway through the splitting procedure, we have numerically 
shown that the fidelity between this output and the output of a fritter into which — 1 
particles are inputted (and none are lost) is 1 when the interactions are ignored {V = 0). 
Even when interactions are included the loss of fidelity is very small and so spontaneous 
emissions will not be the limiting factor in this scheme. 

Taking all these practical limitations into account we find the limit for S is five. All that 
is required to produce these devices is a raising and lowering of the potential barriers and 
so far fewer components are used than in optical schemes to produce devices with the same 
number of ports. 



VI. CONCLUSIONS 



We have proposed a straightforward scheme for implementing atomic beam splitters with 
up to five input and output ports. It requires modulation of the intensity of the optical lattice 
in which the atoms are trapped - something that is readily achievable in the laboratory. 
Importantly, in this scheme, the multiport devices require no more operational complexity 
than an ordinary beam splitter for atoms. 

While experimentally challenging, the prospects for implementing these devices look 
promising and certainly seem to be within reach of current technologies. Indeed lattices 
of the required geometry have already been experimentally demonstrated [l2!]. We have 
considered the effects of potential limiting factors and shown that, in general, they do not 
make the scheme prohibitively difficult. The time the barriers are lowered for requires control 
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to an accuracy that increases with S. However, we found the precisions required to produce 
outputs of the required fidehty could be achieved using current laboratory equipment. The 
effects of several other practical limitations were also found to be minimal. Unfortunately, 
interactions were found to limit the number of ports to five. However, the scheme is simple 
and has the great advantage that, for many applications, the lattice sites do not have to be 
addressed individually, which considerably simplifies its implementation. 

This scheme is also versatile. The inverse transforms, for example, can be achieved 
using the same set-up just by altering the timings. In a similar way, it is possible to achieve 
splitters with unbalanced outputs or variable 'reflectivities'. It should also be straightforward 
to combine multiport devices into useful schemes such as multipath interferometers, and the 
'ring' lattice means they have great potential for use in gyroscopes that can measure angular 
momentum with Heisenberg-limited precision. The versatility of these devices means that 
they are likely to have intriguing prospects for creating interesting entangled states and play 
a useful role in a range of new quantum technologies. 
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FIG. 1: Ring configuration of tlie sites in the optical lattice for the tritter. The lines denote 
coupling between sites due to tunnelling. 




FIG. 2: These plots show how the vary with time for multiport devices with S = 3,4,5 and 6 
from top left to bottom right respectively. We see that for 5 = 3 the Q intersect at t = 2tt/9J and 
for 5 = 4 they intersect at t = tt/4J. For 5 = 5 and 6 there is no intersection of all the 0, in the 
time range, t = — > 27r/J, shown here. 
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Jt/7I 

FIG. 3: Using x ^ 10~^ as our defining criterion for when a balanced splitter is achieved we find a 
five site splitter is produced at t « 5.27r/J. The near intersection of all the O at this time can be 
seen here. 




FIG. 4: The fraction of the total number of atoms, A'^i, N2 and N-^, in each of the three output ports 
of a three-path interferometer given by Eq. (jl9p . The applied phase, (p, inside the interferometer 
is the same across each pair of adjacent sites. 
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FIG. 5: The Jt required to produce a balanced splitter for different S. We find these Jt scale with 
S as 1.8x10-457-0. 
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FIG. 6: The fidelity of the output states of a fritter with interactions (calculated from the overlap 
with the output when there are no interactions) plotted over the range VN/ J = ^ 1 for different 
numbers of atoms. 
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FIG. 7: The correlation between VN/J at a critical fidelity ol F = 0.95 and the number of atoms. 
The crosses are numerically calculated data points for values of up to = 40. The solid curve 
is a line of best fit intended to find the scaling. The scaling is, (1^A^/J)f=o.95 ~ 0.85N^'^'^'^ . 
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FIG. 8: Plot of F/J at a critical fidelity of 0.95 as a function of 5 for = 5. We see V/ J decreases 
rapidly as S increases. For S = 7 the interaction strength required to achieve a fidelity of at least 
0.95 is y/J ~ 10~^. Such a value is likely to be experimentally challenging and so we take 5 = 5 
as the maximum number of sites allowed when interactions are accounted for. 
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FIG. 9: The relationship between the absolute time error at a critical fidelity of F = 0.95 and 
the number of atoms, N, for a fritter. The solid curve is a line of best fit and has the form 
eir=o.95 ~ 0.16N~^'^^ . We see that as the number of atoms increases, the timing error that can be 
tolerated decreases. 
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FIG. 10: Plot of e/r for 5 = 3, 4, 5, 7 and 9 when we input five atoms into the system. We determine 
e/r by taking e at the critical fidelity of 0.95 and dividing it by the Jt required to produce balanced 
splitters (r) for the different S. The fractional time error, e/r, that can be tolerated decreases 
rapidly with S. 
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